Lack of consensus in social systems 
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We propose an exactly solvable model for the dynamics of voters in a two-party system. The 
opinion formation process is modeled on a random network of agents. The dynamical nature of 
interpersonal relations is also reflected in the model, as the connections in the network evolve with 
the dynamics of the voters. In the infinite time limit, an exact solution predicts the emergence 
of consensus, for arbitrary initial conditions. However, before consensus is reached, two different 
metastable states can persist for exponentially long times. One state reflects a perfect balancing 
of opinions, the other reflects a completely static situation. An estimate of the associated lifetimes 
suggests that lack of consensus is typical for large systems. 



Concepts and tools of modern noncquilibrium statis- 
tical physics lend themselves very directly to describing 
complex interacting systems, including phenomena which 
rely on human behavior, e.g. the emergence of collective 
organization in social systems. Recently variants of the 
voter model [3, H 0] have been used intensively to study 
collective phenomena, such as opinion formation or con- 
sensus creation 0, 1, @, S i [I EI EI El EI EI El- 
Many of these efforts have focused on regular lattices 
, which is justified in physical situations, 
but not in the context of the social sciences. In socio- 
cultural situations, the interaction patterns between in- 
dividuals find a better characterization as complex net- 
works in which the connections or relationships (links) 
between individuals (nodes) can change in time. More 
precisely, the full dynamics of such a social network con- 
sists of (i) the opinion formation process taking place on 
the nodes, and (ii) the evolution of the underlying topo- 
logical structure (links). The coupling between these two 
processes reflects how the connections of people influ- 
ence their opinions, and how their opinions determine, 
in turn, their new connections. Although there is an in- 
creasin g re cent interest in mode ling voter dynamics on 
graphs H 11 1 , and on networks 12 , II 14 1 , the dynam- 
ics of the links (ii) is still disregarded in these studies. 
The coevolution of nodes and links - i.e., of the full net- 
work structure - has been studied only in a few works 
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In this letter, we investigate the voter dynamics on an 
adaptive disordered network characterized by (i)-(ii). In 
our network, each node j ("individual") carries a spin 
Oj ("opinion") which can take two different values <r,j — 
±1 [l[. At each time step, (i) the spins are updated 
random sequentially based on a simple majority rule: if 
they are connected to more positive than negative spins, 
their state will be positive in the next time step, and 
negative otherwise; in the case of a tie, the spin remains 
unchanged. Further, (ii) the links are updated as follows: 
two nodes carrying equal (unequal) spins are connected 
with probability p (q). In this letter, we focus on the 
special case q = 1 — p, leaving the general case to [IH . 

As an interpretation, we propose that this model mim- 



ics a two-party electoral system. During a campaign, the 
supporters of one party are keen to interact with sup- 
porters of the other party to try to change their opin- 
ion. This situation can be described by this model with 
p < q, when each agent has more interactions with op- 
ponents than with agents sharing the same opinion (ac- 
cording to the motto that "convinced people do not need 
to be convinced again"). On the other hand, when p > q, 
the agents tend to interact more with individuals sharing 
the same opinion (according to the motto "united we are 
stronger"). The latter behavior seems to be a simplified 
description of the process of political polarization, when 
all the members of a party agree with the official position 
of the party, as often occurs in post-election periods. 

Our model is aimed to describe a free public debate 
in the sense that it does not consider the effects of cen- 
tral institutions or the mass media; neither lobbying, nor 
organized strategies (apart from possibly influencing the 
probability p) are taken into account. As a result, the 
model may also be appropriate to describe groups defined 
by criteria such as education, religion or ethnicity, rather 
than political opinion. Cultural assimilation, the spread- 
ing of a language or a religion of an ethnic or religious 
minority, and social reforms are examples of phenomena 
which can be modeled in this fashion. 

To describe the dynamics of the system, let us focus 
on p(t), the "popularity" of + opinions, defined as the 
average fraction of + nodes at time t. Thus, we consider 
P(M,t), the probability of finding the network with M 
positive spins at time t, and 
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where N is the total (fixed) number of nodes. Clearly, 
p = or 1 correspond to a complete ordering of the 
system, while p = 0.5 characterizes the completely dis- 
ordered state. Contrary to the voter model on regular 
lattices, the global magnetization (m = 2p — 1) is not 
conserved here, but the dynamics is still Z2 symmet- 
ric (i.e., invariant under the global inversion ai t— > — Cj, 
M/N h-> 1-M/N). 
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Since the spins on the nodes flip one at a time, the 
time evolution of P(M, t) is a simple birth-death process 
for which we can write a master equation: 



d t P(M, t) = b M -iP(M - 1, t) + d M+ iP(M + 1, t) 
- [b M + d M ]P(M,t). 



(2) 



Here, &m denotes the birth rate of a positive spin (i.e., 
the rate for flipping a negative to a positive spin), and 
<1m its death rate (flip rate from positive to negative). 
Both depend on M, the current number of positive spins 
in the system. Whether a positive spin will flip or not is 
determined by the number of positive and negative spins 
it is connected to. Since these connections are estab- 
lished randomly, with probabilities p and q, respectively, 
they are controlled by binomial distributions. For exam- 
ple, the probability that a positive spin is connected to 
exactly k of the other M — 1 positive spins, is given by 
B M -i, P (k) = ( M ^ 1 )p k (l-p) M - 1 - k . Writing a similar 
expression for the probability of this spin to be connected 
to k negative spins, the death rate is given by 
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The prefactor simply reflects the probability to find a 
positive spin among the N spins. The step function 
Q(k' — k) expresses the fact that if the selected spin is 
connected to k positive and k! negative spins then, in the 
next time step, it will take the state sgn(fc — k'). Simi- 
larly, the birth rate is: 
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For the special case q = 1 — p [2l| , the 0-functions can 
be eliminated from bu and dm, due to the properties of 
the binomial distributions, to yield: 
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Before solving the master equation, we discuss its ap- 
proximate solutions in the thermodynamic limit, in order 
to identify the different types of behavior and the param- 
eter regimes where they occur. For N — > oo, the bino- 
mial distribution B>N yP (k) approaches a normal distribu- 
tion with mean Np, so that cLm is given by the Gaussian 
error function multiplied by the prefactor M/N , 
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0, if M > Np, 
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FIG. 1: The death dm and birth dm rates (continuous and 
dashed lines, respectively), for N = 1000 and different prob- 
abilities p. For fixed p, there are three different regimes of 
behavior bounded by Eq. ©. 



apart from a region of width yj Np(l — p) around Np. 
Similarly, the birth rate &m is described by the comple- 
mentary error function. These forms of the transition 
rates (see Fig. [IJ, along with the probability p and the 
initial fraction of positive spins, Mq/N, determine the 
late-time properties of the model. The master equation 
controls the flow of M/N, as a function of time, leading 
to four distinct regimes, depending on the relative mag- 
nitudes of Mq/N and p. In a {p, M /N) phase diagram, 
these different regimes are bounded by 
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For p < 0.5 and M /N < p, we find that M/N stays 
below p at later time also. Indeed, in the approximation 
([5]), we have a pure death process 

d t P(M, t) = (M + 1) P{M + 1, t) - MP (M, t) , (7) 

which leads to the extinction of the positive population. 
The steady state, = \im.t^nx p(t) — 0, is reached ex- 
ponentially as p{t) ~ po exp(— i/iV). Similarly, if p < 0.5 
and Mq/N > 1 — p, we have a pure birth process, and 
the system relaxes exponentially to the state p^, — 1 on 
the same characteristic time scale as in the previous case 
(due to the Z2-symmetry). In the intermediate region 
Mq /N e \p, 1 — p] , the dynamics is described by: 



Nd t P(M,t) = (M + 1)P(M + 1, t) 
+[N - M + 1)P(M - 1, t) 



NP(M,t), (8) 



and the system reaches a disordered phase: p^ = 0.5. 
Again, the relaxation is exponential, with a characteristic 
time scale N/2. 

For p > 0.5 the pure death and pure birth regimes are 
the same as for p < 0.5. A small minority (Mo/iV < 
1 — p) will become extinct, a large majority iMo/N > p) 
will win. However, a new feature appears in the interval 
Mq/N G [1 — p,p], where the system seems to acquire 
infinite memory. Both the death and birth rates vanish 
in this region, so that dtP(M, t) — whence the fraction 
of positive spins remains frozen at its initial value, Mq /N. 

Our analytic findings are tested by simulations 2J| on 
a network with N — 1, 000 nodes. The relaxation into the 
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FIG. 2: Time evolution of the fraction M/N of positive spins 
for a network of N = 1000 nodes, for p = 0.3 (a) and p = 
0.7 (b). The values are averaged over 1000 runs. Time is 
measured in Monte Carlo steps, tM-C ~ ispm-flip/A^ . The insets 
show the exponential relaxation to the final states. 




FIG. 3: The fraction of positive spins M/N after tMC = 10, as 
function of p and Mo/N for a network of N = 1, 000 nodes in 
3D (a), and 2D (b) displaying clearly the phase transitions. 



four late-time states is displayed in Figs. [2ja) for p < 0.5 
and (b) for p > 0.5. The possible outcomes of the voter 
dynamics, for all parameters p and initial fractions Mo /N 
of the positive population are summarized in Fig. [3] 

To illustrate the picture further, Fig. [5] shows the 
outcome of the voter dynamics for two initial fractions 
of positive population: one starting from a minority 
Mq/N < 0.5, and one starting from a majority Mo/N > 
0.5. For small p, the system reaches a disordered state, 
independent of Mq/N: the "open mindedness" of the 
population (reflected by a large probability 1 — p to com- 
municate with the opposite party) leads to an equal dis- 
tribution of opinions. In contrast, an "inflexible attitude" 
(characterized by a large probability p of linking up with 
similar opinions) leads to an unchanging distribution of 
opinions. For intermediate values of p, the system reaches 
a completely ordered state: all voters reach the same 
opinion, namely that of the initial majority. 
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FIG. 4: Outcome of the voter dynamics as function of p. Cross 
sections of Fig. taken along M /N = 0.3 and M /N = 0.7. 



In conclusion, in the thermodynamic limit, the voter 
dynamics has four possible outcomes: a perfect balance 
of opinions, a static situation, or consensus (p = 0, 1). In 
this last section, we discuss how these findings are modi- 
fied in finite systems. The two completely ordered states 
are absorbing states, thus they will be reached from the 
other two (metastable) states. Two interesting questions 
remain: First, in which of the two absorbing states will 
each metastable state arrive, and second, how do the re- 
laxation times depend on system size? 

To answer the first question, we write Eq. ^ in a 
matrix representation, dt \v(t)) = L \v(t)), where \v(t)) is 
the (N + l)-dimensional column vector with components 
P(M, t), M = 0, 1, N, and the time evolution operator 
L can be read off from Eq. @ . The steady states of the 
system are the eigenvectors of L with zero eigenvalues: 
|0 ) = (1, 0, 0) T and |0jy) = (0, 0, l) T . To 
find their adjoints, it is convenient to study the symmet- 
ric/antisymmetric states: |0±) = (|0rj) ± |0jv))/2. Im- 
posing orthonormality (0±|0±) = S±, we obtain from 
(0±|L = the right eigenvectors (0+| = (l, 1, l) 
and (0_| = (l, aci, — £2, — x i, — l), with 

_ ((( r i+i + l Vj+2 + i)r j+ 3 + ••• + l)r„_i + 1 
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Thus we can compute explicitly the 
final state |V>oo) = liiD-t— Kx^t)- Given an initial state 
|V>o) } the solution to Eq. ©is \tp t ) = J2» e_Af '*l/ i )(A t l^o) J 
where and \fj) are the left and right eigenvectors of L 
corresponding to eigenvalues A M . Expecting no other zero 
A M 's, we have |Voo) = |0+)(0+|Vo) + |0_)(0_|Vo). For 
example, if the initial state is a population with Mo < 
N/2 positive spins, the final state is 



~ (l+x Mo , 0, 



0, l-x 



(9) 



This result shows that, indeed, in the t — > 00 limit, the 
system arrives in one of its two absorbing states. More- 
over, it provides the relative probabilities with which ei- 
ther will be reached, as a function of the initial Mq. Fur- 
ther details will be published elsewhere 2l| . 
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FIG. 5: The probability P(N, 00) = (1 - x M )/2, to reach the 
state with M — N positive spins at infinite times: P(N, 00) = 
implies certain extinction of the positive population, while 
P(N, 00) — 1 represents a purely positive population. 

These findings are confirmed in Fig. [51 obtained by 
direct iteration of the master equation. For p > 0.5, the 
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FIG. 6: Time r needed to reach consensus as function of the 
system size for p = 0.3, M /TV = 0.4 (a), and p = 0.8, M /TV = 
0.4 (b). The insets show the exponential dependence on TV. 



dynamics reduces to a simple majority rule: The final 
state is completely ordered, following the opinion of the 
initial majority. For p < 0.5, a small minority disappears, 
while a big majority wins the competition, just as in the 
case of infinite system size. A novel effect occurs for 
intermediate values of the initial positive population. In 
finite systems, the fully disordered (metastable) state has 
a finite lifetime, during which the system forgets its initial 
condition. Then, after a very long time, it randomly falls 
into one of its absorbing states. To rephrase, the initial 
positive population is equally likely to become extinct or 
to take over the whole system. In this regime, the final 
outcome of the voter dynamics is completely random. 

Finally, we explore numerically the relaxation times 
into the absorbing states in Fig. [51 The lifetime r of 
the metastable states increases with the system size TV, 
as t ~ e a (p) Jv j with a p-dependent coefficient a(p). An 
estimate of the time to consensus for a network of TV = 
1000 voters with p = 0.3 shows that it may take as many 
as 10 36 spin flips to reach one of the absorbing states. In 
much larger systems, consensus is practically impossible. 

In conclusion, the best strategy for a minority group 
is to establish many contacts with its opponents. In this 
way, it can convince half of them and keep this balance for 
a long time. If the same group is less open for discussions, 
it cannot overcome the majority, but at least it will not 
disappear. It is tempting to speculate how these results 
might be applied to real social systems. Will two-party 
systems, once formed, persist for long times? Will bilin- 
gual regions remain bilingual? Will relatively isolated 
parties continue to receive the same, almost constant per- 
centage of the vote? Will closed religious communities 
continue to exist without gaining or losing members? 

These particular results are obtained from a simple 
adaptive model in which certain important social factors 
are neglected (e.g., spatial and age structures, a spectrum 
of opinions, etc.). More precise statements will certainly 
need additional assumptions on the character of the net- 
work or opinion dynamics. Nevertheless, even this simple 
model can give a better understanding of empirical data 
not yet explained. For example, a study of the num- 
ber of languages in the Solomon Islands [23j found that 
small islands (less than 100 square miles) have a single 



language, but above this size the number of languages 
increases. The finite size effect pointed out by our model 
can be a possible explanation of this phenomenon. 

The importance of the model presented here stems 
from the fact that it is mathematically simple, exactly 
solvable, and easily generalized to more complex situ- 
ations (for q ^ 1 — p see [2l|, for epidemics networks 
see [IS] ). Our model provides a method how to describe 
analytically the adaptive nature of the interpersonal re- 
lations, and by this, it can serve as a "baseline model" 
which captures the key characteristics of social systems, 
namely, having disordered networks of agents and dy- 
namically changing connections between them. 
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